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RAPPORT  V' EVALUATION  TECHNIQUE 

VU  SYMPOSIUM  ORGANISE  PAR  LA  COMMISSION  V E  OVNAM1QUE  VES  FLU  WES 
OE  L‘ AGARD  SUR 

LE  CALCUL  VE  L' INTERACTION  FLUIVE  PAREAIT-ELUIVE  V1SQUEUX 


J.C.  LE  BALLEUR 

0(iEce  National  d'Etudu  it  de  Recfie'tc/iw  A(A06patialei  (ONE RA) 
29,  avenue  de  la  V iviiion  Leclinc  -  9232 0  CHATILLON  -  FRANCE 


I.  INTROPUCTION  - 

Du  29  septembre  au  ler  octobre  1980,  la  Commission  de  Dynamlque  des  Fluides  de  l’ACARD,  prdsldde  par  le 
Dr.  Orlik-Ruckemann,  a  organise  un  Symposium  sur  le  "Calcul  de  1 ’Interaction  Fluide  Parfait-^luide  Visqueux" 
dans  les  locaux  de  l'Alr  Force  Academy  des  Etats-Unis  3  Colorado-Springs,  Colorado,  USA.  La  grande  qualite 
de  sa  preparation  technique  et  matjrielle  a  dtl  ressentie  par  tous  les  participants. 

Le  Symposium  Stait  organise  par  un  Comlte  de  Programme  international,  dirige  par  M.  l'Ing.  en  Chef 
B.  MONNFRIE  et  par  le  Dr.  B.  QlIINN  qui  ont  egalement  contribud  3  la  prdsldence  des  Sessions  et  de  la  Discus¬ 
sion  finale. 

Les  quatre  conferences  gdnerales  et  les  27  communications  presentees  sont  representatives  des  recherches 
numdrlques  orlentees  vers  le  calcul  des  ecoulements  de  fluide  visqueux  3  grands  nombres  de  Reynolds  au  moyen 
de  methodes  dltes  "interactives",  reallsant  un  couplage  entre  le  calcul  des  couches  vlsqueuses  et  le  calcul 
du  fluide  parfait  externe. 

Les  communications  ont  en  pratique  ete  regroupdes  dans  les  trois  sessions  sulvantes  : 

I  -  Ecoulements  non  decolles  -  Concept  de  couche  mince. 

II  -  Forte  interaction  en  regime  turbulent  sans  decollements  etendus. 

III  -  Ecoulements  decolles. 

Comme  l'ont  montrd  le  detail  des  exposes  et  la  discussion  de  Table  Ronde,  ces  trois  types  de  preoccupa¬ 
tions  ont  en  fait  ete  presentes  dans  de  nombreuses  communications  au  cours  des  trois  sessions.  Cette  reparti¬ 
tion  constltue  done  davantage  une  indication  sur  trois  aspects  d'un  meme  probleme  encore  lmparfaltement 
rdsolu,  notamment  en  ce  qui  conceme  les  questions  de  forte  interaction  et  de  ddcollement,  plutot  que  la 
presentation  de  trois  voles  de  recherche  distinctes. 

L'ensemble  des  textes  associds  aux  Conferences  et  aux  communications  est  public  dans  l'"AGARD  Conference 
Proceeding  n"  291",  "Computation  of  Viscous-Inviscid  Interactions".  Ce  document  dtant  deja  disponible  pour 
une  dtude  approfondle  des  exposes,  alnsi  que  de  la  discussion  de  Table  Ronde,  l'objectif  du  present  rapport 
se  limits  3  une  analyse  globale  succincte  des  communications  presentees,  ainsi  qu'3  l'examencfes  conclusions 
et  des  recommendations  qui  peuvent  s'en  degager. 


2.  ANALYSE  GENERALE  - 

2.1  -  NiviaiLK  d' appAoxAJncuU.cn  du  a alculA  - 

Les  motivations  de  ce  Symposium  etaient  d'dvaluer  les  possibilit6s  qui  seront  offertes  3  court  et  moyen  termes 
par  les  methodes  numeriques  pour  le  calcul  et  1 'optimisation  des  performances  adrodynamiques,  dans  l'hypothese 
d'une  analyse  des  dcoulements  en  fluide  visqueux  ne  s'appuyant  pas  sur  une  resolution  directe  des  equations 
de  Navier-Str kes  moyenndes  appliquees  a  tout  le  fluide.  Ce  souci  correspond  3  l'idde  que  le  cout  des  resolu¬ 
tions  globales  "Navier-Stokes"  sur  des  maillages  adaptds  aux  grands  nombres  de  Reynolds, et  1 ' lncert itude  des 
modules  de  turbulence  mis  en  jeu, feront  de  ces  mdthodes  un  outil  de  prevision  opdrationnel  3  plus  long  terme, 
notamment  en  dcoulement  trid imensionnel. 

"ne  evaluation  succincte  des  approximations  numeriques  utilisables  en  fluide  visqueux  laminaire  ou  turbulent 
aux  grands  nombres  de  Reynolds  est  donnde  par  la  Conference  du  Pr.  Kline  [22]  et  par  celle  d ' introduct ion 
gendrale,  Le  Balleur  [1].  Leur  synthese  conduit  en  resume  3  degager  les  4  -iveaux  d 'approximation  suivants  : 

1  -  Fluide  parfait 

2a  -  FLuide  Parfait  +  correlations  emplrlques  adfmensionnelles 
2b  -  Fluide  Parfait  +  couches  limites  (couplage  faible) 

2c  -  Approximations  parabolisees  des  equations  de  Navier-Stokes 
2d  -  Fluide  Parfait  +  zones  decollees  isobares 

3a  -  Fluide  Parfait  +  couches  visqueuses  minces  (couplage  fort) 

3b  -  Fluide  Parfait  +  zones  Navier-Stokes  (couplage  fort) 

l.  -  Resolutions  directes  globales  des  equations  de  Navier-Stokes  : 

4a  -  Equations  movennees  +  module  de  turbulence 

4b  -  Equations  filtrdes  (simulation  des  grosses  structures  turbulentesl . 

4c  -  Equations  completes  (simulation  numertque  de  la  turbulence). 


II  faut  ajouter  d  cea  approximations  num^riques  les  Etudes  analytlques  compiymentai res,  asymptot lques  ou 
non-asymptotiques,  dgvelopp^es  pour  lea  problfcmes  d ' interact  ion  couche  liraite-onde  de  choc  ou  de  bord  de 
fuite,  qui  peuvent  etre  rattacWes  au  niveau  d 'approximat Ion  3  des  mythodes  de  couplage  fort. 


La  plupart  des  approximations  pratiquement  utilises  reinvent  actuellement  des  nlveaux  1  et  2«  Le  b^n^f Ice 
des  progr£s  r^cents  realises  dans  les  mythodes  nvun^riques  relatives  aux  ecoulements  de  fluide  parfait  pourralt 
etre  utilise  au  niveau  3,  qui  repr^sente  le  seuil  h  partir  duquel  tous  les  ph^nomdnes  de  forte  interaction 
visqueuse  contenus  dans  les  Equations  de  Navier-Stokes  moyennees  peuvent  etre  mod£lis£s.  au  moins  qualita- 
t ivement . 

Le  niveau  3a  correspond  aux  object  if s  du  present  Symposium.  Sa  mtse  S  1’etat  op^rationnel  est  attendue  par 
le  groupe  de  Stanford  pour  le  courant  de  la  prochaine  d^cennie  [221.  Ce  niveau  3a  lnclut,  conune  les  nlveaux 
3b  et  4,  la  possibility  de  calculer  les  ecoulements  avec  ou  sans  dycollements  dans  la  limite  de  precision 
des  equations  de  couche  mince  adoptees,  la  model Isat ion  allant  des  Equations  de  Prandtl  aux  Equations  de 
Navier-Stokes,  objets  des  methodes  de  couplage  du  niveau  3b.  Le  Symposium  montre  que  1 'etat  de  1 'art  se 
situe  actuellement  3  un  stade  intermed ia Ire  entre  les  nlveaux  2b,  2d,  3a. 

2.2  -  AppioxAmatA-OMi  intz>im$dicU'iei>  zntAe  couplage  ficubte  zt  couplage  jo>tf  (niveau  2)  - 

Compte  tenu  de  1 ’orientation  du  symposium  vers  les  methodes  numyriques,  peu  de  presentations  ont  £te 
consacryes  au  niveau  2a  et  3  Sexploitation  de  correlations  expyr imentales ,  mentionnyes  cependant  par 
Cersten  et  al  [23]  pour  les  problemes  de  recollements  subsonlques.  Tmplicitement ,  les  corryiations  emplri- 
ques  sont  egaleraent  presentes  dans  les  model i sat  ions  les  plus  simples  dyveloppyes  pour  1  *  interact  ion  onde  de 
choc-couche  limite, en  regime  turbulent  transsonique,  dans  lesquelles  on  admet  que  le  fluide  parfait  est 
perturbe  par  une  "rampe  visqueuse  yquivalente"  dont  1 ’angle  de  dyflexlon  maximale  se  situe  sur  le  plan  de 
l’hodograohe  entre  la  deviation  maximale  et  la  deviation  sonique.  Les  progr£s  numeriques  effectuys  sur  cette 
approche  ont  ety  indiquys  par  Wai,  Yoshihara  C 9 ] ^  Jou ,  Murman  [15].  Une  approximation  de  meme  nature  est 
utilisee  par  Desopper,  Crenon  [5].  File  est  citye  par  Stanewsky  et  al.  [4]  pour  1 ’exploitat ion  de  la  methode 
de  Inger  [IS].  Implicitement ,  le  recours  &  des  correlations  exper imentales  est  encore  prysent  dans  certaines 
relations  de  fermeture  utilisees  pour  le  calcul  des  couches  limites, ainsi  que  dans  les  modeles  simples 
de  nlveaux  2b,  2d,  qui  sont  frequemment  incorpores  pour  restituer  les  effets  globaux  des  petlts  dycol lements 
de  bord  de  fuite,  ou  des  bulbes  laminaires  courts  de  bord  d'attaque  [6,  19,25]. 

Aux  nlveaux  d 'approximat ion  2b,  2d,  les  couches  visqueuses  sont  calculees  a  partir  des  yquations  de  Prandtl, 
resolues  pour  une  distribution  de  pression  externe  fixee.  Elies  sont  couplees  a  1 ’ycoulement  de  fluide 
parfait  externe  par  des  procedures  approchees,  dites  de  couplage  faible.  Ces  procedures  maintiennent  une 
influence  dominante  du  fluide  parfait  externe  sur  la  solution  finale,  soit  en  vertu  de  theories  asymptot i- 
ques  de  faible  interaction,  soit  en  raison  de  1  ’emploi  d ’algor ithmes  de  couplage  fort  non  converges,  soft 
encore  par  des  techniques  de  lissage  arbitraires  sur  l'effet  de  deplacement.  Les  principaux  obstacles  de  ces 
approches  sont  la  gyneration  bien  connue  de  solutions  singulieres  au  voisinage  du  decollement,  alnsi  que  la 
suppression  du  phenomene  d’ influence  amont  en  ecoulement  visqueux  supersonique ,  c*est-a-dire  en  pratique 
1' inaptitude  au  calcul  et  au  positionnement  des  interactions  couche  limite-onde  de  choc  fl). 

Sur  les  singularites  de  decollement  et  leur  suppression  par  des  methodes  de  resolution  inverses, dans  le  cas 
des  methodes  integrales  d’entralnement,  une  svnthese  des  resultats  analvtiaues  acauis  par  les  deux  eroupes 
de  1*0NERA,[1]  [14],  sur  les  equations  bidimensionnelles  stat ionnaires  ou  instationnaires,  et  tr idimension- 
nelles  stat Ionnaires,  est  donnee  par  Cousteix,  Houdeville  [14],  Elle  rappelle  que  1 'ut il isat ion  de  methodes 
de  couche  limite  a  pression  externe  imposee  comme  indicateur  de  decollement  ajoute  aux  approximations  du 
modele  de  calcul  les  insuf f isances  d'une  resolution  mathemat iquement  hasardeuse.  II  semble  cependant  que 
cette  difficulty  puisse  parfois  etre  surmontee  dans  1 ’approximat ion  de  niveau  2d,  qui  correspond  a  la  mody- 
lisation  non-visoueuse  de  decollement s  isobares  pour  les  regimes  de  decrochaee.  ainsi  ou’en  tymolgnent  les 
resultats  de  Carlson  [26]  pour  les  profils  et  voilures  subsoniques.  De  meme,  en  ecoulement  bid imensionnel 
instat ionnaii e,  Desopper  et  Crenon  [5]  demontrent  la  possibility  d'obtenir  des  solutions  de  couplage  faible 
jusqu'au  voisinage  du  decollement  sur  des  profils  ou  des  pales  d ’helicopteres  en  ecoulement  incompressible  ou 
transsonique.  En  ecoulement  tr id imensionnel  stationnaire  sur  des  corps  elances,  Fiddes  [30]  demontre  par 
contre,  en  etudiant  les  ecoulements  coniques,  que  la  determination  des  1 ignes  de  decollement  a  l'origine  des 
nappes  tourbillonnaires  au  moyen  de  methodes  de  couche  limite  traditionnelles  de  faible  interaction  n'est 
probablement  plus  acceptable,  meme  si  les  lignes  de  syparation  sont  reconsiderees  de  maniere  iterative. 
L’addition  d’une  methode  de  forte  interaction,  couche  limite  couplee  ou  presentement  [30]  modele  asymptot ique 
en  triple  couche,  semble  etre  indispensable. 

Dans  le  cas  mieux  connu  des  ycoulements  bid imensionnels  stat ionnaires,  qui  ont  servi  de  support  a  la  quasi- 
totalite  des  presentations  de  ce  Symposium,  1’ importance  du  recours  a  des  methodes  de  couplage  fort  pour  la 
determination  correcte  des  points  de  decollement  est  illustree  par  les  compara isons  de  calculs  de  couche 
limite  "interact if s"et  de  calculs  de  couche  limite  classique,  comme  celles  effectuees  par  Gersten  et 
al.  [23]  pour  1 'ycoulement  incompressible  sur  une  marche  descendante  emoussee,  ou  par  Fortunato  [27] 
dans  le  cas  d'une  rampe  de  compression  supersonique. 

Le  souci  d'effectuer  des  calculs  de  couche  limite  interactifs  est  evident  dans  les  mythodes  de  calcul 
classiques  des  profils  d’ailes,  ou  les  techniques  de  resolution  de  1 'equation  du  potentiel  ineluent  une 
iteration  directe  sur  un  calcul  d'ypaisseur  de  deplacement.  Bien  que  ces  mythodes  interdisent  une  evaluation 
non-empirique  des  dycollements,  et  que  les  techniques  numeriques  de  couplage  effectivement  mises  en  jeu  ne 
garantissent  pas  une  capacite  systymatique  3  resoudre  les  phynomenes  de  couplage  fort,  de  bord  de  fuite  ou 
de  sillage,  divers  progr^s  dans  cette  approche  sont  mis  en  evidence  par  le  Symposium.  Le  dyveloppement  des 
methodes  de  singularities  est  utilise  par  Oskam  [191,  Butter  et  Williams  [25]  pour  le  calcul  des  profils 
multi-corps  en  incompressible.  L' extension  des  methodes  de  relaxation  pour  1  Equation  complete  du  potentiel 
est  indlquee  par  Rosch  et  Klevenhusen  [7],  par  Leicher  [6],  pour  le  calcul  de  profils  multi-corps  en  regime 
transsonique.  Une  technioue  inverse  de  genyration  de  profils  supercr it iques  sans  choc  avec  effet  de  couche 
limite  est  donnee  par  Nebeck  et  al .  [31. 
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Enfin,  Stanewsky  et  al.  [4]  utillsent  la  method®  de  Inger  [181  pour  amillorer  la  provision  da  l*(palaaaur 
de  deplacement  dans  1 1  Interact  Ion  couche  llmlte-onde  de  choc  sur  profile  transsonlques,  la  couplage  au 
fluide  parfalt  etant  toutefols  liss£  par  une  representation  polynomlale  du  deplacement. 


Dans  la  llmlte  de  validity  des  equations  retenues  pour  le  calcul  de  I'epalsseur  de  deplacement,  la  conver¬ 
gence  d'une  iteration  sur  I'epalsseur  de  deplacement  peut  conduire  &  la  realisation  d'une  methode  de 
couplage  fort  du  niveau  d 'approximat ion  3a.  II  suffit  essent iellement  [1],  que  d'une  part  la  relation  de 
couplage  soit  consistante  et  depourvue  de  module  de  llssage,  et  que  d 'autre  part  lea  echelles  de  discreti¬ 
sation  numerique  solent  sufflsamment  fines,  en  cas  de  bulbes  de  decol lenient  par  exemple. 

Une  premiere  conclusion  importante  pour  la  classification  des  methodes  de  calcul  du  fluide  vlsqueux  conslste 
done  H  dissocier  le  concept  de  couplage  faible  et  celui  d'effet  de  deplacement.  Les  conferences  generales 
du  Symposium,  Lock  [2],  Le  Balieur  [1],  attlrent  de  surcrolt  1 'attention  sur  la  necesslte  de  ne  pas  asslmller 
le  concept  de  deplacement  de  forte  interaction  h  l'idde  taclte  d 'approximat ion  de  couche  mince,  d'hypoth&ses 
de  couche  llmlte,  ou  d'equations  de  Prandtl.  11  ressort  en  particulier  [1]  [2]  que,  moyennant  le  calcul 
d'un  ecoulement  de  fluide  parfalt  fictif  superpose  a  l'ecoulement  reel  au  sein  des  couches  vlsqueuses,  11 
est  possible  de  definir  rigoureusement  le  deplacement  ainsi  que  la  relation  de  couplage  jusqu'au  niveau 
d 'analyse  3b,  correspondent  au  couplage  des  equations  de  Navier-Stokes. 

Dans  les  methodes  de  couplage  fort  du  niveau  3a,  les  approximations  de  couche  mince  qui  sont  retenues  pour 
les  equations  des  couches  vlsqueuses  se  repercutent  pour  une  part  sur  la  precision  du  calcul  de  l'effet  de 
deplacement.  Elies  alterent  de  plus  1* evaluation  du  champ  de  presslon  dans  les  couches  vlsqueuses  et  par 
consequent  le  couplage  en  effet  de  courbure  des  sillages  [2]  [10]  [1],  notamment  au  volslnage  des  bords  de 
fuite.  Les  progres  en  ces  domaines  sont  recherches  pour  une  part  dans  la  confrontation  £  1 'experience  de 
modeles  mathemat iques  faisant  appel  tantot  3  une  couche  limite  unique  [1]  [2]  [11]  [13],  tantot  S  une 
ou  plusieurs  sous-couches  [10]  [17]  [18].  Les  progrds  sont  recherches  pour  une  autre  part  dans  la  determi¬ 
nation  de  developpements  asyraptot iques  rationnels  pour  les  solutions  des  equations  de  Navier-Stokes  &  la 
limite  infinie  du  nombre  de  Reynolds,  dans  l'hypothese  d'une  extrapolation  possible  aux  nombres  de  Reynolds 
usuels  [10]  [16]. 

Pour  les  modeles  utilisant  une  couche  limite  unique,  la  conference  d ' introduct ion  gen£rale  [1]  developpe 
la  possibilite  de  traiter  les  gradients  de  pression  normaux  des  couches  vlsqueuses  en  s'appuyant  sur  une 
formulation  deficitaire  des  equations  par  rapport  a  un  ecoulement  de  fluide  parfalt  fictif  superpose  . 

Le  orocede  elimine  en  particulier  les  anomalies  de  couches  supercritiques  au  sens  de  Crocco-Lees,  meme 
si  l'on  se  limite  S  un  niveau  d 'approximat ion  du  premier  ordre.  II  apparaft  par  allleurs  souhaitable  [1] 

[2]  que  les  corrections  vlsqueuses  de  "courbure",  qui  donnent  acces  aux  approximations  d'ordre  superieur 
du  champ  de  presslon,  solent  evaluees  5  partlr  de  la  courbure  moyenne  des  lignes  de  courant  indulte  dans 
le  fluide  parfalt  par  l'effet  de  deplacement  de  forte  interaction.  Un  calcul  couple  du  proche  slllage 
vlsqueux  respectant  sa  dissymgtrie  semble  de  ce  fait  indispensable  l 7].  Cette  model isat ion  const itue  le 
fondement  theorique  des  analyses  de  couplage  fort  utilisant  des  methodes  integrates  pour  calculer  les  couches 
vlsqueuses  [1]  [11]  .  Au  premier  ordre  d 'approximation ,  une  technique  de  meme  nature  est  indiqu^e  par 
Wilmoth  et  Dash  [13]  pour  evaluer  les  effets  vlsqueux  des  probl&mes  d 'arr iere-corps  avec  jets. 

Sur  les  modeles  asymptot iques  multi-couches  developpes  pour  les  probl&mc;,  oe  oord  de  fuite  ou  d '  interact  Ion 
couche  limite  -  onde  de  choc  en  regime  transsonique  turbulent,  une  excellente  revue  S  caractdre  synthetique 
est  donnee  par  la  conference  generate  de  Melnik  [10]  ,  que  complete  la  presentation  de  Adamson  et 
Messiter  [16].  Supposant  que  la  theorle  asymptot ique  de  couplage  faible  reste  valable  en  dehors  des  domaines 
de  forte  interaction  dont  elle  determine  les  conditions  inltiales,  les  modeles  asymptot iques  locaux  corapie- 
mentaires  lndiquent,  sauf  rare  exception,  1' importance  du  gradient  de  pression  normal  dans  la  couche 
visqueuse  turbulente,  ainsi  que  le  role  relativement  passif  des  sous-couches  vlsqueuses.  Les  resultats 
analytlques  de  Melnik  et  al.  pour  les  bords  de  fuite  minces  et  exempts  de  decollement  peuvent  etre  couples 
£  une  methode  de  couche  limite  Interactive  dans  les  calculs  de  profils  trans9onlques  [10].  En  revanche, 
l' exploitation  des  modeles  asymptotiques  pour  1 'Interact Ion  sous  le  choc  reste  plus  difficile,  dans  le  cas 
des  profils  transsoniques ,  dans  la  mesure  ou  les  effets  vlsqueux  reduisent  le  niveau  global  de  recompression 
et  elolgnent  de  la  situation  asymptotique  de  choc  fort. 

Les  modules  non-asymptot iques  a  deux  couches  se  rattachant  aux  concepts  deLighthill  sont  £voqu£s  par  la 
Conference  de  Melnik  [10]  ainsi  que  par  les  presentations  de  Bohning  et  Zierep  [17]  et  de  Inger  [18].  De 
meme  que  dans  les  theories  asymptotiques  il  s'aglt  de  modules  locaux,  utilisables  seulement  au  volslnage  des 
interactions  couche  limite  -  onde  de  choc,  et  negligeant  la  viscosite  dans  la  partie  externe  de  la  couche 
visqueuse.  Une  difference  essentielle  reside  par  contre  dans  le  role  determinant  accorde  icl  a  la  sous- 
couche  visqueuse.  Les  resultats  recents  [17]  [18]  sont  donnes  sous  forme  d'etude  parametrlque  Incluant  les 
problemes  de  courbure  de  paroi  et  de  decollement  naissant. 

2.4  -  CalcuZb  Incluant  dtcoltejnent*  de  couche  mince.  ~ 

Les  modules  multi-couches  dlsponibles  pour  les  dcoulements  turbulents,  asymptotiques  ou  non-asymptotiques, 
perdent  leur  validlte  en  cas  de  decollements,  meme  s'ils  se  llmltent  3  de  petlts  bulbes.  Presentement , 
seuls  les  modules  a  couche  limite  unique  permettent  d'aborder  le  calcul  des  decollements  turbulents. 

La  modellsation  la  plus  simple  repose  sur  la  resolution  "interactive"  des  equations  de  Prandtl  et  d'un 
ecoulement  potent iel  externe  raccorde.  Elle  sera  it  consistante  dans  le  cas  laminaire  avec  le  module  asympto¬ 
tique  en  triple  couche.  Des  resultats  de  ce  type  sont  recherches  pour  des  ph^nomenes  d *  interact  ion  locale,  en 
ecoulement  incompressible  [121  l 20]  [21]  ou  supersonique  d'onde  simple  [27]  [28]  [29],  la  perturbation  de 
l'ecoulement  externe  de  l'ecoulement  pouvant  alors  etre  calculge  d'une  faqon  explicite,  ou  par  1  *  intermedia ire 
d* integrates  de  Cauchv-  Cebeci  et  al.  (20]  etudient  le  decollement  au  bord  d'attaque  d'un  profil  elliptique 
en  incidence  et  obtiennent  des  solutions  en  analogic  avec  l'edatement  des  petits  bulbes  laminaires  ou 
transltlonnels.  Gleyzes  et  al.  [21]  s'appuient  sur  une  etude  exper imentale  detailiee  de  bulbes  transit ionnels 
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au  bord  d'attaque  d'un  profil  pour  definir  un  module  empirique  approche,  permettant  un  calcul  de  la 
transition  dans  les  bulbes  courts  avec  effet  de  couplage  local.  Ardouceau  et  al .  [28]  calculent  le 
decollement  turbulent  sur  une  rampe  supersonlque  par  une  technique  iterative  utlllsant  une  resolution 
Inverse  des  equations  de  Prandtl.  La  discussion  num£rique  du  choix  optimal  de  la  methods  Inverse  conflrme 
1' importance  du  report  il  la  parol  de  la  condition  de  couplage,  procede  connu  pour  eiiminer  les  comporte- 
ments  supercritiques  111. 

Des  modules  de  forte  interaction  3  couche  limite  unique  plus  g£n£raux,  capables  d'un  calcul  des  regions  de 

decollement,  sont  en  cours  de  d^veloppement  pour  les  ecoulements  complexes  autour  de  profils  II)  19]  (11). 

Des  calculs  transsonlques  dans  1 'approximation  des  petites  perturbations  sont  indlqu^s  par  Wai  et 
Yoshlhara  [91,  le  decollement  au  pied  du  choc  restant  toutefois  assujetti  a  une  modelisation  de  rampe 
visqueuse  equivalente.  Des  calculs  fondes  sur  la  resolution  de  1' equation  complete  du  potent  lei, coupiee 
S  une  methode  lntegrale  qul  modelise  l'dcart  entre  fluide  parfait  fictif  et  flulde  reel, ont  pu  etre 
realises,  Le  Balleur  [“].  Tls  incluent  le  calcul  et  le  positlonnement  des  slllages  dissyraetr iques ,  le  calcul 
des  bulbes  de  decollement  sur  le  profil  ou  au  bord  de  fuite,  3  l'exception  des  decollements  sous  le  choc  en 
regime  transsonique.  Les  progres  obtenus  par  Neron  dans  la  definition  d'une  methode  de  slngularlt^s  numerique- 
ment  bien  condit ionnee,dans  le  cas  des  profils  a  bord  de  fuite  anguleux  ou  ef files,  ont  permis  d* adapter  la 
methode  precedente  au  calcul  de  profils  multi-corps  avec  bulbes  de  decollements  multiples,  en  ecoulement 
Incompressible  111]. 

2.5  -  Calcul  de6  pKo^it. 6  el  de*  cuXah  -  M£thodc6  num&fuqueA  de  couplage  - 

Les  methodes  de  calcul  des  profils  transsonlques  ayant  atteint  un  stade  de  developpement  plus  operationnel 

prennent  en  compte  glneralement  les  effets  de  deplacement  et  de  courbure  du  slllage,  mais  ne  sont  pas  en 
mesure  de  calculer  les  decollements.  Des  elements  d'analyse  sont  indiques  dans  les  conferences  de  Melnik 

[10]  ,  Le  Balleur  M]«  Une  revue  plus  complete  a  ete  donnee  par  la  Conference  de  Lock  [2],  presentee  par 
le  Dr  Green.  Cette  revue  couvre  egalement  les  methodes,  nettement  moins  nombreuses,  qul  ont  ete  developpees 
pour  les  voilures  transsonlques,  sous  forme  de  methodes  iteratives  entre  des  calculs  potentiels  et  des 
calculs  de  couche  limite  tridimensionnels.  En  raison  des  difficultes  restant  a  resoudre  en  ce  domaine  il 
n'est  pas  encore  possible  de  deboucher  sur  des  methodes  de  couplage  fort.  Neanmolns,  ainsi  que  le  demontre 
la  presentation  de  Firmin  [8]  sur  des  configurations  voilure-fuselage  en  theorie  des  petites  perturbations 
transsonlques,  une  prise  en  compte  tres  appreciable  de  la  viscosite  peut  deja  etre  realisee. 

Les  methodes  de  calcul  operationnelles  precedentes,  destinies  auxprofils  et  voilures  transsonlques  en 
I'absence  de  decollement,  s'appuient  pour  realiser  le  couplage  numerique  entre  fluide  parfait  et  couches 
visqueuses  sur  un  calcul  iteratif  direct  de  1 'effet  de  deplacement,  dont  la  convergence  est  d'autant  moins 
assuree  que  le  decollement  est  plus  voisin  ou  que  le  maillage  de  calcul  est  plus  fin.  Dans  ces  methodes,  la 
stability  du  calcul  est  au  mieux  obtenue  par  des  techniques  de  sous-relaxat ion  relat ivement  empiriques. 

Les  progres  recents  invitent  au  contraire  a  considerer  comme  un  Element  central  des  methodes  de  couplage  fort 
la  definition  mathematique  d 'algorithmes  de  couplage  numeriquement  stables.  Une  analyse  du  probleme,  ainsi 
qu'une  revue  des  principales  methodes  numeriques  deja  developpees  pour  le  couplage  fort,  sont  donnees  dans 
la  Conference  d' introduction  generale,  Le  Balleur  HI*  Dans  le  cas  bidimensionnel  stationnaire ,  les  methodes 
de  relaxation  directes  peuvent  etre  stabilisees  sans  empirisme  en  calculant  un  coefficient  de  sous-relaxat ion 
local,  tandis  que  des  techniques  de  relaxation  semi-inverses  ou  semi-implicites  permettent  d'aborder  le 
couplage  fort  des  couches  limites  decollees.  Les  presentations  de  Wai  et  Yoshihara  [91/  Le  Balleur  et  Neron 

[11] ,  Veldman  [12],  Ardonceau  et  al.  [28]  completent  1 'analyse  du  probleme. 

2.6  -  Aot*ea  It/pe*  d'appsiocheA  - 

A  l'Scart  des  autres  methodes  de  calcul  evoquees  dans  ce  Symposium,  la  presentation  de  Ashurt  et  al,  [24] 
foumit  une  tentative  devaluation  de  methodes  lagrangiennes  utlllsant  des  particules  tourbillonnaires 
discretes  pour  simuler  numeriquement  les  decollements  turbulents  de  grande  dimension.  Appliquee  au 
recollement  derriere  une  marche  descendante  et  restreinte  a  des  hypotheses  de  calcul  bidimensionnelles,  bien 
qu'une  a t tray ante  Simula t ion  qualitative  soit  obtenue,  cette  technique  ne  semble  pas  en  mesure  de  recouper 
l'experience  pour  la  determination  des  tensions  de  Reynolds,  ou  plus  simplement  pour  reproduire  la  variation 
de  la  longueur  de  recollement  en  fonction  du  nombre  de  Reynolds. 

Il  semble  enf in  interessant,  comme  peut  le  sugg^rer  la  conference  ggnerale  de  Kline  [22],  de  replacer  le 
deroulement  du  present  Symposium  dans  le  cadre  plus  vaste  du  calcul  des  ecoulements  turbulents  complexes, 
qui  a  ete  choisi  pour  theme  de  la  conference  de  Stanford  1980-81.  Celle-ci  se  propose  de  determiner  une 
gamme  d 'exper imentat ions  detaillees  couvrant  des  ecoulements  turbulents  aussi  divers  que  possible,  et  de 
les  utiliser  pour  eprouver  les  methodes  numeriques  actuellement  disponibles. 


3  -  C0MMENTA7RES  PART7CU17ERS  - 

Des  remarques  plus  ponctuelles  peuvent  etre  degagees  de  ce  Symposium  en  s'appuyant  sut  l'exemple  du  calcul 
des  profils  d'ailes,  qui  a  constitue  pour  une  bonne  part  le  centre  d'interet  des  presentations.  Cet  exemple 
cumule  en  effet  la  plupart  des  problemes  d 1  in teract ion  visqueuse  de  1 'aerodynamique  externe  transsonique 
bidimensionnelle,  et  constitue  un  prelude  aux  developpements  de  calculs  tridimensionnels. 


3.  /  -  UtltiA&tion  de.  V Equation  du  potdntczl  - 

La  totalite  des  presentations  s'appuie  pour  determiner  1 ‘ Seoul ement  de  fluide  parfait  externe  sur  la  resolu¬ 
tion  de  liquation  du  potentiel,  eventuellement  simplifiee  dans  le  cadre  de  1 'approxlmat ion  des  petites 
perturbations  transsonlques.  Des  progres  sont  done  encore  possibles  en  ce  qui  concerne  le  calcul  de 
1'pooulement  bidimensionnel  du  fluide  parfait.  Cette  remarque  est  d£j3  confirmee  par  un  certain  nombre  de 


tentatives  pr£sent£es  dans  ce  Symposium  pour  am£liorer  le  cond i t ionnement  numgrique  des  methodes  de  singula¬ 
rity,  utilises  pour  les  profils  multiples  aux  basses  '’itesses.  l>a  remarque  s’applique  en  outre  essentiel- 
lement  3  1 ’absence  d *ut il lsat ion  des  mgthodes  num£riques  r£solvant  les  Equations  d'P.uler,  probablement  en 
raison  de  leur  cout  plus  elev£,  bien  que  les  nombres  de  Mach  locaux  mis  en  Jeu  en  flulde  vlsqueux  transsonlque 
au  voislnage  du  decollement  naissant  sous  le  choc  rendent  d§j3  contestable  1 ’approximat ion  potentlelle. 

Cette  premiere  limitation  d’origine  non-visqueuse  dolt  etre  presente  3  1 ’esprit  dans  les  comparisons  3 
1 ’experience.  File  a  par  allleurs  donne  naissance  3  plusleurs  ecoles  de  pens£e  pour  1 'ut 11 lsat Ion  pratique 
des  methodes  numeriques  potentielles  transsoniques,  fondees  sur  Sexploitation  de  techniques  num£riques 
tantot  conservatives,  tantot  non-conservatives,  tantot  par t iel lement  conservatives  vis-3-vis  du  debit 
masse.  Cette  confusion  est  d^autant  plus  delicate  qu'elle  interfere  g£n£ralement  aver  1 ’adoption  d'un 
traitement  plus  ou  moins  elabore  pour  le  couplage  des  effets  vlsqueux  au  pied  du  choc,  le  r£sultat  global 
relevant  alors  d’une  certatne  compensation  d’erreur  entre  les  deux  effets.  Cet  obstacle  theorique  s’applique 
surtout  aux  techniques  non-conservatives,  reputees  pour  leurs  resultats  plus  volsins  de  1 'experience  en 
1 'absence  de  tout  effet  vlsqueux.  Les  techniques  conservatives  convenablement  coupl£es  aux  effets  vlsqueux 
sous  le  choc  paraissent  a  ce  jour  offrir  une  vole  plus  sure  pour  1 'evaluation  difficile  de  la  trainee. 


Utitt^ation  dz  m$thodz.A  u\t^.gKaZzA  pout  Zz  zaZcuZ  dz- 4  couche  ZajhZZ&a 


t’ne  seconde  remarque  se  degageant  des  calculs  de  profils  ou  de  voilure  est  le  recours  generalise  3  des 
methodes  integrates,  le  plus  souvent  a  des  methodes  d 'ent ralnement ,  pour  calculer  les  couches  visqueuses 
turbulentes.  Si  cette  remarque  ne  s'applique  pas  a  toutes  les  analyses,  notamment  a  certaines  etudes  de 
phenomenes  vlsqueux  locaux,  si  elle  peut  s’expliquer  en  premier  lieu  par  le  souc i  de  reduire  le  cout  du 
calcul  numerique,  notamment  en  vue  d'une  extension  aux  problemes  tr id imensionnels  ,  il  est  clair  cependant 
que  ce  cholx  se  justifie  par  la  quallte  des  resultats  pratiques  qui  peuvent  etre  obtenus  pour  les  profils 
de  vitesse  moyenne  et  pour  le  deplacement.  Cette  chance  de  simplification  numerique  paratt  en  outre  pouvoir 
conduire  au  calcul  des  bulbes  de  decollement ,  ainsi  qu'a  1 ' incorporation  d'un  calcul  d'hysteresis  sur  les 
contraintes  turbulentes  pour  ameliorer  la  modelisation  de  1 'entrafnement .  He  plus,  les  progres  contenus  dans 
une  formulation  deficitaire  des  equations  visqueuses  par  rapport  a  un  ecoulement  de  fluide  parfait  fictif 
superpose, (1 ],  mettent  en  evidence  le  fait  que  les  methodes  integrales  constituent  en  realite  un  modele  de 
calcul  moins  restrictif  que  les  equations  de  Prandtl.  Leur  utilisation  couplee  avec  la  resolution  d’un 
ecoulement  de  fluide  parfait  fictif  prolongeant  1 'ecoulement  externe  conduit  en  effet  a  une  representation 
approchee  des  gradients  de  pression  normaux  plus  conforme  aux  equations  de  Navier-Stokes  des  couches 
minces.  Elle  tolere  en  particulier  l’apparition  d'ondes  de  choc  internes  aux  couches  limites. 


iz  COUAbuAZ  dZA  AA 


L* influence  de  1 'effet  de  displacement  des  sillages  sur  le  calcul  de  1 'ecoulement  global  est  unanimement 
admise.  Par  contre,  le  consensus  apparent  sur  1 ’ importance  de  l'effet  de  courbure  des  sillages  masque 
en  realite  des  formulations  differentes,  en  fonction  de  la  technique  adoptee  pour  definir  une  courbure 
movenne  des  lignes  de  courant.  Si  les  diverses  formulations  sont  sensiblement  eauivalentes  pour  le  sillage 
lointain  qui  n'exerce  qu'une  faible  influence,  un  traitement  d lssymetr ique  du  proche  sillage  est  par 
contre  probablement  necessaire  a  la  realisation  correcte  d'un  calcul  de  forte  interaction,  ce  calcul  suppo- 
sant  en  outre  l'emploi  d'un  algorithme  de  couplage  stable  et  depourvu  de  lissage.  Ces  difficultes  expliquent 
vraisemblablement  qu'un  veritable  consensus  ne  soit  pas  encore  atteint  pour  la  prise  en  compte  des  effets 
visqueux  de  bord  de  fuite  en  regime  turbulent,  a  l'exception  du  passage  a  la  limite  infinie  du  nombre  de 
Reynolds.  Les  consequences  sur  1 'evaluation  pratique  de  la  trainee  semblent  pouvoir  etre  importantes. 


3.4  -  IntzAaction  choz-couckz  ZatmXz 

Les  avis  divergent  aussi  sur  le  choix  d'un  traitement  pratique  de  1  *  interact  ion  couche  limite  -  onde  de  choc 
adapte  au  calcul  des  profils  ou  des  voilures  transsoniques .  Aux  incertitudes  dela  mentionnees  sur  la  technique 
numerique  adoptee  pour  I'ecoulement  externe,  s'aioute  la  difficult^  d'un  maillage  de  discretisation  generale- 
ment  trop  grossier  pour  le  phenomene,  dont  l'etendue  totale  est  souvent  inferieure  a  la  maille  de  calcul. 

Bien  qu'il  soit  generalement  admis  que  les  consequences  globales  de  ce  phenomene  sur  I'ecoulement  soient  moins 
determinantes  que  celles  resultant  de  la  region  de  bord  de  fuite,  les  imperatifs  des  calculs  sont  au  moins 
de  deux  natures.  Le  calcul  visqueux  approche  doit  d'une  part  modifier  fondamentalement  1 'evolution  de  la 
recompression  sous  le  choc  et  en  reduire  le  niveau,  ce  oui  exclut  1 'utilisation  direete  d'un  modele  asvmpto- 
tique  avec  choc  droit  et  qui  peut  exiger  dans  certains  cas  l'addition  d'une  technique  de  rampe  visqueuse 
art  if ic ielle.  Le  calcul  approche  doit  d 'autre  part  estimer  1 'evolution  globale  de  la  couche  limite  avec  une 
precision  suffisante  pour  calculer  valablement  la  region  de  bord  de  fuite.  Le  developpement  futur  de 
resolutions  a  echelles  plus  fines,  garantlssant  un  calcul  et  un  posit ionnement  du  phenomene,  paralt  toutefois 
indispensable,  notamment  en  cas  de  decollement  local. 
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En  depit  de  ces  difficultes  non-resolues,  les  methodes  de  calcul  des  profils  transsoniques,  fondees  sur  la 
resolution  de  1 'equation  du  potentiel  couplee  avec  des  methodes  integrales  pour  les  couches  visqueuses,  ont 
atteint  un  niveau  operationnel  de  prevision  pratique,  au  moins  en  1 'absence  de  decollement.  Les  performances 
obtenues  sont  superieures  a  celles  des  methodes  de  resolution  direete  des  equations  de  Navier-Stokes  moyennees 
a  la  fots  en  cout,  en  quallte  et  en  densite  de  maillage.  Ces  resultats  encourageants  en  ecoulement 
bidimens ionnel  stationnaire  ne  doivent  cependant  pas  masquer  les  lacunes  evoquees  plus  haut ,  ainsi  que  la 
difficulte  qu'il  y  a  a  ce  niveau  d 'approximat ion  pour  effertuer  des  comparaisons  objectives  entre  les 
calculs  et  1 ' experimentation  transsonlque.  11  serait  en  particulier  souhaitahle  de  disposer  d ' exper i ences 
de  references  assez  purement  bid imens ionnel les,  c 'est -3-d ire  3  peu  pres  exemptes  «  'effets  de  parols, 
comportant  un  sondage  tres  precis  des  couches  limites  en  aval  du  ddclenehement  de  la  transition,  dans 
lesquelles  la  comparalson  au  calcul  pourrait  etre  effectuee  pour  un  nombre  de  Mach  et  une  incidence  determine 
sans  ambiguite,  plutot  que  par  1 'artifice  tres  contestable^ d 'un  ajustement  de  la  portance  globale. 
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Les  progr^s  recents  effectues  sur  les  modeles  mathemat iques  couplant  fluide  parfait  et  couches  visqueuses, 
ainsl  que  sur  les  algorithmes  numeriques  permettant  de  resoudre  Its  p robl ones  de  couplage  fort  correspondant s , 
ont  attelnt  le  stade  permettant  une  extension  des  methodes  bid Imensionnel les  stat ionnalres  prec£dentes  aux 
regimes  des  decollements  de  couches  minces.  Dans  ce  domaine,  qui  correspond  aux  grands  coefficients  de 
portance,  aux  prof  11s  multi-corps  ou  au  buffeting  transson ique,  une  sensibillte  accrue  aux  model isat ions 
turbulentes  alnsi  qu'au  traitement  de  la  zone  de  transition  doit  toutefois  etre  attendue  dans  les  perfor¬ 
mances  pratiques  des  methodes  de  calcul.  Par  ailleurs,  un  lien  theorique  entre  ces  modeles  de  calcul  du 
decollement  et  du  recollement,  et  les  modeles  plus  simples  a  lignes  de  }et  isobares  des  grands  decollements 
paratt  devoir  etre  recherchee  dans  le  futur. 

Bien  que  l'essentiel  de  l1 etude  reste  3  accomplir,  l’extension  des  methodes  de  couplage  fort  s'appuyant  sur 
une  resolution  en  fluide  parfait  de  l'ecoulement  externe  et  sur  des  methodes  integrales  pour  les  couches 
visqueuses  devrait  pouvoir  deboucher  sur  le  calcul  pratique  des  ecoulements  tr id imensionnel s  ou  instat ionna 1  res 
avec  decollements  de  couche  mince.  A  cet  effet,  1 'obstacle  theorique  et  numerique  le  plus  nouveau  a  franchir 
consiste  sans  doute  a  imaginer  un  modele  visqueux  de  forte  interaction  decrivant  le  lachcr  et  1 'enrouletaent 
des  nappes  de  sillage  tourbillonna ires  issues  des  corps  fuseles. 

4  -  CONCLUSIONS  ET  REC0MMANVAT10NS  - 

Le  Symposium  demontre  clairement  que  1* effort  consent!  ces  dernieres  annees  dans  le  domaine  des  methodes 
numerlques  de  calcul  du  fluide  parfait  peut  etre  largemenl  valorise  d'un  point  de  vue  pratique  par  1 ‘addition 
d'une  resolution  complements! re  interactive  des  couches  visqueuses. 

Le  travail  est  largement  avance  en  ce  qul  concerne  les  ecoulements  bidimensionnels  stationnaires  utilisant  une 
approximation  potentielle  de  l'ecoulement  externe.  L’etat  de  1 'art  se  situe  en  ce  domaine  a  un  niveau  inter- 
mediaire  entre  des  methodes  de  couplage  faible,  des  methodes  de  modelisation  isobare  pour  les  grands  decolle¬ 
ments,  et  des  methodes  de  couplage  fort,  ces  dernieres  <?tant  capables  de  model iser  1 'ensemble  des  phenomenes 
reels  d ' interaction  visqueuse,  au  moins  d'une  maniere  qualitative. 

Les  tendances  des  progres  observes  vont  dans  le  sens  des  methodes  de  couplage  fort,  avec  le  developpement 
conjoint  de  modeles  mathemat iques  couplant  fluide  parfait  et  couches  visqueuses,  ainsi  que  d 'algorithmes 
numerlques  indispensables  a  leur  resolution  couplee  rigoureuse.  Bien  que  des  incertitudes  demeurent  sur  les 
problemes  de  bord  de  fuite  ou  d *  interact  ion  couche  limite-onde  de  choc,  des  methodes  operat ionnell es  et 
performantes  existent  pour  les  profils  transsoniques  ;  on  peut  observer  notamment  1 'apparition  de  methodes 
nouvelles  capables  de  calculer  les  decollements  de  couches  minces  ou  de  bord  de  fuite. 

La  mise  en  oeuvre  de  resolutions  des  equations  d 'Euler  completes  pour  1 'eeoulement  externe  paraitrait  souhai- 
table,  d'une  part  pour  eliminer  le  cboix  delicat  entre  les  techniques  numeriques  potentielles  conservatives 
ou  non-conservatives,  d'autre  part  pour  etendre  le  domaine  d 'appl icat ion  des  methodes  de  couplage  fort  aux 
problemes  d 'aerodynamique  interne  transsonique. 

Le  travail  est  beaucoup  moins  developpe  dans  le  domaine  des  ecoulements  instat ionnalres  ou  tr id imensionnel s . 
Neanmoins  les  premiers  resultats  obtenus  par  resolution  de  1 ' eeoulement  de  fluide  parfait  dans  1 'a- proximal i on 
potentielle  et  par  addition  d'un  calcul  couple  des  couches  visqueuses  au  moyen  de  methodes  integrales 
paraissent  demontrer  la  possibility  de  generaliser  les  modeles  bidimensionnels. 

Les  performances  part icul i erement  encourageantes  des  methodes  de  couplage  fort  en  vue  des  applications 
pratiques,  pour  lesquelles  elles  constituent  un  outil  de  simulation  numerique  approchee  des  equations  de 
Navier-Stokes  en  couches  minces,  conduit  a  recommander  pour  le  futur  d'associer  plus  systema t iquement  et 
plus  etroitement  le  developpement  de  ces  methodes  numerlques  et  celui  des  techniques  do  resolution  des 
equations  du  fluide  parfait. 

De  meme,  les  developpements  futurs  et  complementa ires  des  methodes  de  couplage  fort  ainsi  que  Jes  techniques 
globales  de  resolution  des  equations  de  Navier-Stokes  auraient  probabl ement  interet  a  etre  plus  etroitement 
assoc  les.  II  en  est  ainsi  par  exemple  sur  les  problemes  concernant  la  limitation  des  dornaines  de  calctil 
"Navier-Stokes"  aux  seules  Regions  visqueuses.  II  apparalt  en  outre  possible  de  dTvelopper  de  nouveaux 
algorithmes  pour  les  equations  de  Navier-Stokes,  au  moins  dans  1 'approximat ion  des  couches  minces,  en 
s'appuyant  sur  les  techniques  numeriques  issues  des  methodes  de  couplage  fort. 
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I.  INTRODUCTION  - 

From  29th  September  to  1st  October  1980,  the  AGARD  Fluid  Dynamics  Panel,  chaired  by  Dr  K.  Orlick  Ruckemann,  organized 
a  Symposium  on  the  “Computation  of  Viscous-lnviecid  Interaction*  at  the  US  Air  Force  Academy.  Colorado  Springs,  Col.  US.A 
The  high  quality  of  its  technical  and  material  preparation  was  appreciated  by  all  participants. 

The  symposium  was  organized  by  an  international  Program  Committee,  headed  by  Chief  Engineer  8.  MONNERIE  and  Dr.  8.  QUINN, 
who  also  chaired  sessions  and  the  final  discussions. 

The  four  general  lectures  and  27  papers  presented  surveyed  current  numerical  research  oriented  towards  the  calculation  of 
viscous  fluid  flows  at  high  Reynolds  numbers  by  means  of  so-called  "interactive"  methods,  realizing  a  coupling  between  the  com¬ 
putation  of  viscous  layers  and  that  of  the  external  ideal  fluid. 

The  papers  were  in  fact  grouped  into  the  three  following  sessions  : 

I  -  Unseparated  flows  -  Thin  layer  concept. 

II  -  Turbulent  strong  interaction,  without  extensive  separated  flow  regions. 

Ill- Separated  flows. 

As  was  shown  in  the  detailed  papers  as  well  as  in  the  Round  Table  Discussion  preoccupations  concerning  these  three  categories 
were  evident  in  many  papers  during  each  of  the  three  sessions.  This  categorization  thus  indicates  three  aspects  of  the  same  pro¬ 
blem  still  imperfectly  solved  (especially  as  regards  strong  interaction  and  separation)  rather  than  three  distincts  avenues  of 
research. 

The  complete  set  of  texts  associated  with  the  Lectures  and  Papers  is  published  in  the  AGARD  Conference  Proceedings  CP  291, 
"Computation  of  Viscous-lnviscid  Interaction".  This  document  being  already  available  for  a  thorough  study  of  the  papers,  as 
well  as  of  the  Round  Table  discussion,  the  purpose  of  this  Report  is  limited  to  a  brief  overall  analysis  of  the  papers  presented, 
as  well  as  a  survey  of  the  conclusions  and  recommendations  that  can  be  deduced  therefrom. 

3.  GENERAL  ANALYSIS  - 

3.1  -  Approximation  levels  of  computations 

The  purpose  of  this  symposium  was  to  evaluate  the  possibilities  offered  at  short  and  medium  term  by  numerical  methods  for  the 
computation  and  optimization  of  the  aerodynamic  performances  of  aircraft  within  the  assumption  of  an  analysis  of  viscous 
flows  that  is  not  based  on  a  direct  solution  of  the  averaged  Navier-Stokes  equations  applied  to  the  whole  flow  field.  This  viewpoint 
corresponds  to  the  idea  that  the  cost  of  overall  "Navier-Stokes"  solutions  on  mesh  patterns  adapted  to  high  Reynolds  numbers, 
and  the  uncertainty  of  the  turbulence  models  involved,  will  make  these  methods  an  operational  prediction  tool  at  a  longer  term, 
in  particular  in  three-dimensional  flow. 

A  brief  evaluation  of  the  numerical  approximations  usable  in  laminar  or  turbulent  viscous  fluid  at  high  Reynolds  numbers  is 
given  by  the  Prof.  Kline  Lecture  [22]  and  by  the  General  Introduction  Lecture,  Le  Balleur  [1].  Their  synthesis  leads  one,  in 
short,  to  bring  to  light  the  four  following  levels  of  approximation  : 

1  -  Ideal  fluid 

2a  -  Ideal  fluid  +  nondimensional  empirical  correlations 
2b  -  Ideal  fluid  +  boundary  layers  Ivwak  coupling) 

2c  -  Parabolized  approximations  of  the  Navier-Stokes  equations 

2d  —  Ideal  fluid  +  isobaric  separated  zones 

3a  —  Ideal  fluid  +  thin  viscous  layers  (strong  coupling) 

3b  —  Ideal  fluid  +  Navier-Stokes  zones  (strong  coupling) 

4  -  Direct  overall  solutions  of  the  Navier-Stokes  equations  : 

4a  —  Averaged  equations  +  turbulence  model 

4b  —  Filtered  equations  (simulation  of  the  large  turbulent  structures) 

4c  -  Complete  equations  (numerical  simulation  of  turbulence). 

To  these  numerical  approximations  should  be  added  complementary  analytic  studies,  whether  asymptotic  or  not,  developed  for  the  pro¬ 
blems  of  shock-boundary  layer  or  trailing  edge  interaction.  These,  in  turn,  may  be  related  to  the  third  level  of  approximation  of 
strong  coupling  methods. 
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Molt  of  the  approximations  commonly  used  at  present,  concern  levels  1  and  2.  The  advantages  of  the  recent  progrea  achieved 
in  numerical  methods  applied  to  ideal  fluid  flows  could  be  used  at  level  3.  This  would  constitute  a  threshold  from  which  all 
the  phenomena  of  strong  viscous  interaction  contained  in  the  averaged  Navier-Stokes  equations  might  be  modelled,  at  least 
qualitatively. 

Level  3a  corresponds  to  the  objectives  of  the  present  Symposium.  It  is  expected  to  be  brought  into  operational  status  by  the 
Stanford  team  within  the  next  decade  (22] .  This  level  3a,  as  well  as  levels  3b  and  4,  include  the  possibility  of  computing 
flows  without  or  with  separations,  within  the  accuracy  limits  of  the  thin  layer  equations  adopted,  the  modeling  extending 
from  the  Prandtl  equations  to  the  Navier-Stokes  equations,  subjects  of  the  coupling  methods  of  level  3b.  The  Symposium 
shows  that  the  state  of  the  art  lies,  at  present,  at  an  intermediate  stage  between  levels  2b.  2d  and  3a, 

2.2  -  Intermediate  approximations  between  weak  coupling  and  strong  coupling  ( level  21 

Considering  the  orientation  of  the  Symposium  toward  numerical  methods,  few  papers  were  devoted  to  level  2a  and  to  the  exploi¬ 
tation  of  experimental  correlations,  mentioned  however  by  Gersten  et  al.  (23)  for  the  problems  of  subsonic  reattachment  . 
Implicitly,  empirical  correlations  are  also  present  in  the  simplest  models  developed  for  the  shock  wave-boundary  layer  interac¬ 
tions,  in  turbulent  transonic  regime,  in  which  it  is  assumed  that  the  fluid  is  perturbed  by  an  "equivalent  viscous  wedge"  whose 
maximum  deflexion  angle  is  located,  in  the  hodograph  plane,  between  the  maximum  deflexion  and  the  sonic  deflexion.  The  nu¬ 
merical  progress  achieved  in  this  approach  has  been  described  by  Wai  and  Yoshihara  [9]  Jou  and  Murman  [15].  An  approxima¬ 
tion  of  the  same  nature  is  used  by  Desopper  and  Grenon  [5  J.  It  is  mentioned  by  Stanewsky  et  al.  [ 4[  for  the  exploitation  of 
the  Inger  method  [18].  Implicitly,  recourse  to  experimental  correlations  is  still  present  in  some  closure  relationships  used  for  the 
computation  of  boundary  layers,  as  well  as  in  simple  models  of  levels  2b,  2d,  which  are  frequently  incorporated  to  account  for 
the  overall  effects  of  small  trailing  edge  separations,  or  of  the  small  laminar  bubbles  at  the  leading  edge  |6.  19.  25]. 

At  approximation  levels  2b,  2d,  the  viscous  layers  are  computed  from  the  Prandtl  equations,  solved  for  a  given  external  pressure 
distribution.  They  are  matched  to  the  external  ideal  fluid  flow  by  approximate  processes,  of  so-called  weak  coupling.  These  pro¬ 
cesses  maintain  a  dominant  influence  of  the  external  ideal  fluid  on  the  final  solution,  either  because  of  asymptotic  theories  of 
weak  interaction,  or  through  the  use  of  non-converged  strong  coupling  algorithm^  or  again  by  arbitrary  smoothing  techniques  on 
the  displacement  effect.  The  principal  obstacles  to  these  approaches  are  the  well-known  generation  of  singular  solutions  in  the 
vicinity  of  separation,  as  well  as  the  suppression  of  the  upstream  influence  phenomenon  in  viscous  supersonic  flow,  ie.  in  prac¬ 
tice.  their  inaptitude  for  computing  and  positioning  shock  wave-boundary  layer  interactions  [1]. 

On  the  separation  singularities  and  their  suppression  by  inverse  solution  methods,  in  the  case  of  integral  methods  of  entrainment, 
a  synthesis  of  the  results  acquired  by  the  two  ONERA  groups, [1],  [14],  on  two-dimensional  steady  or  unsteady,  and  steady 
three-dimensional  equations,  is  given  by  Cousteix  and  Houdeville  [14],  It  recalls  that  the  use  of  boundary  layers  with  assigned 
external  pressure  as  separation  indicator  adds  to  the  approximations  of  the  computation  model,  the  insufficiencies  of  a  hazardous 
mathematical  solution.  It  seems,  however  that  this  difficulty  might  sometimes  be  surmounted  in  the  approximation  of  level  2d, 
which  corresponds  to  the  inviscid  modeling  of  isobaric  separations  for  stalling  regimes,  as  witnessed  by  the  results  of  Carlson  (26) 
for  transonic  profiles,  and  of  Maskewetal.  [31]  for  subsonic  profiles  and  wings.  In  the  same  manner,  Desopper  and  Grenon  [5] 
prove  the  possibility  of  unsteady  weak  coupling  solutions  up  to  the  vicinity  of  separation  on  airfoils  or  helicopter  blades  in 
incompressible  or  transonic  flow.  In  steady,  three-dimensional  flow  on  slender  bodies,  Fiddes  [30]  shows  on  the  other  hand, 
while  studying  conical  flows,  that  the  determination  of  the  separation  lines  at  the  origin  of  vortex  sheets  by  means  of  ronven 
tional  boundary  layer  methods  with  weak  interaction  is  probably  not  acceptable,  even  if  the  separation  lines  are  reconsidered  in 
an  iterative  way.  The  addition  of  a  strong  interaction  method,  coupled  boundary  layer  or  presently,  (30)  a  triple  deck  asymptotic 
model,  seems  to  be  indispensable. 

In  the  better  known  case  of  steady  two-dimensional  flows,  which  served  as  a  support  to  almost  all  papers  of  this  Symposium, 
the  importance  of  using  strong  coupling  methods  for  correctly  determining  reattachment  points  is  illustrated  by  comparisons  of 
"interactive"  boundary  layer  computations  and  conventional  boundary  layer  computations,  such  as  those  carried  out  by  Gersten 
et  al.  [23]  for  the  incompressible  flow  on  a  rounded  backward-  facing  step,  or  by  Fortunato  [27]  in  the  case  of  a  supersonic 
compression  ramp. 

The  problem  of  performing  interactive  boundary  layer  computations  is  obvious  in  the  conventional  methods  for  computing  wings 
profiles,  where  the  techniques  for  solving  the  potential  equation  include  a  direct  iteration  on  the  boundary  layer  displacement 
thickness.  Although  these  methods  preclude  a  non-empirical  evaluation  of  the  separations,  and  the  numerical  coupling  techniques 
actually  implemented  do  not  guarantee  a  systematic  capacity  for  solving  the  strong  interactions,  the  trailing  edge  or  the  wake 
phenomena,  some  progress  in  this  approach  is  brought  to  light  by  the  Symposium.  The  development  of  panel  methods 
is  used  by  Oskam  [19],  Butter  and  Williams  [25],  for  the  computation  of  multicomponent  airfoils  in  incompressible  flow.  The 
extention  of  relaxation  methods  for  the  complete  potential  equation  is  indicated  by  Roach  and  Klevenhusen  |7),  by  Leicher  |6], 
for  the  computation  of  multicomponent  airfoils  in  transonic  regime.  An  inverse  technique  for  generating  supercritical  profiles 
without  shock,  with  boundary  layer  effect,  it  given  by  Nebeck  et  al.  (3).  Lastly,  Stanewsky  et  al-  [4]  use  the  Inger  method 
1 18)  to  improve  the  prediction  of  the  displacement  thickness  in  the  shock  wave  -  boundary  layer  interaction  on  transonic  air¬ 
foils,  the  coupling  to  the  ideal  fluid  being,  however,  smoothed  by  a  polynomial  representation  of  the  displacement. 
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Search  for  strong  coupling  approximations  (level  31 


Within  the  validity  limits  of  the  equations  retained  for  the  computation  of  the  displacement  thickness,  the  convergence  of  an 
iteration  on  displacement  thickness  may  lead  to  the  realization  of  a  strong  coupling  method  of  approximation  level  3a  Essentia 
ly,  1 1 ) ,  it  is  sufficient  that  on  the  one  hand,  the  coupling  relationship  be  consistent  and  deprived  of  smoothing  model,  and  on 
the  other  hand,  the  numerical  discretization  scales  be  sufficiently  line,  in  case  of  separation  bubbles,  for  instance 

A  first  important  conclusion  for  the  classification  of  viscous  flows  computation  methods  thus  consists  of  dissociating  the  weak  coupling 
concept  from  that  of  displacement  effect  The  general  lectures  of  the  Symposium,  Lock  (21,  Le  Balleur  (1|,  emphasize  moreover  the  neces¬ 
sity  of  not  assimilating  the  strong  interaction  displacement  concept  with  a  tacit  approximation  of  thin  layer,  of  boundary  layer  assump 
lions,  or  of  Prandtl  equations.  It  appears  in  particular,  ( 1,  2  (,  that  provided  the  computation  of  a  fictitious  ideal  fluid  flow  is  superposed  to 
the  real  flow  within  the  viscous  layers  it  is  possible  to  define  strictly  the  displacement  as  well  as  the  coupling  relationship  up  to  the 
analysis  level  3b,  corresponding  to  the  coupling  of  the  Navier-Stokes  equations. 

In  the  strong  coupling  methods  of  level  3a.  the  thin  layer  approximations  that  are  retained  for  the  viscous  layer  equations  have  a  first  impact 
on  the  precision  of  the  displacement  effect  computation  Moreover,  they  alter  the  evaluation  of  the  pressure  field  in  the  viscous  layers  and 
consequently  the  coupling  "curvature  effect"  of  the  wakes  (2|,  (101  [U  particularly  in  the  vicinity  of  the  trailing  edges  Progress  in  this  field 
is  sought  in  part,  by  the  comparison  with  experiments  ofmathematical  models  calling  upon  either  a  single  boundary  layer  IH |2L 11 U  |13)or 
one  or  several  sublayers  (10j  (1(|18L 


Progress  is  also  sought  by  the  determination  of  rational  asymptotic  development  for  the  solutions  or 
the  Navier-Stokes  equations  in  the  limiting  case  of  an  infinite  Reynolds  number,  within  the  assumption  of  a  possible  extrapolation  to 
usual  Reynolds numbers(10)[16].  For  the  models  using  a  single  boundary  layer,  the  introductory  general  lecture  |1  (develops  the  possibility 
of  treating  the  normal  pressure  gradient  of  the  viscous  layers,  on  the  basis  of  a  "defect  formulation"  of  the  equations  relative  to  a  superposed 
fictitious  flow  of  ideal  fluid.  The  process  eliminates  in  particular,  the  anomalies  of  supercritical  layers  in  the  sense  of  Crocco- Lees,  even  if  it  is 
limited  to  an  approximation  level  of  the  first  order.  It  also  appears  desirable  [1]  (21  that  the  viscous  correctionsof  "curvature",  which  give  access 
to  the  higher  order  approximations  of  the  pressure  field,  be  evaluated  from  the  mean  curvature  of  the  streamlines  induced  in  the  ideal  fluid  by 
the  strong  interaction  displacement  effect.  Acouplad  computation  of  the  viscous  near  wake  respecting  its  dissymmetry  hence  seems  indispen 
sable  (1 J.  This  modeling  constitutes  the  theoretical  foundation  of  strong  coupling  analyses  using  integral  methods  for  computing  viscous  layers 
(11(11].  At  the  first  order  of  approximation,  a  technique  of  the  same  nature  is  mentioned  by  Wilmoth  and  Dash  ( 13)  for  evaluating  the  viscous 
effects  of  afterbodies  with  jet. 


For  the  asymptotic  multilayer  models  developed  for  the  problems  of  trailing  edge  or  shock  wave-boundary 
layer  interaction  in  turbulent  transonic  regime,  an  excellent  survey  of  a  synthetic  character  isgiven  by  the  general  lecture  of  Melnik  |10],  com¬ 
plemented  by  the  Adamson  and  Messiter  paper  (16|.  Assuming  that  the  asymptotic  theory  of  weak  coupling  remains  valid  outside  the  domains 
of  strong  interaction  of  which  it  determines  the  initial  conditions,  the  complementary  local  asymptotic  models  show,  except  for  rare  excep 
tions,  the  importance  of  the  normal  pressure  gradient  within  the  turbulent  viscous  layer,  as  well  as  the  relatively  passive  role  of  the  viscous  sub 
layers.  The  analytical  results  of  Melnik  et  al.  for  thin  trailing  edges  exempt  of  separation  may  be  coupled  with  a  method  of  interactive  bounda 
ry  layer  in  the  computation  of  transonic  airfoils  (10(.  On  the  contrary,  the  exploitation  of  asymptotic  models  for  the  interaction  under 
shock  remains  more  difficult,  in  the  case  of  transonic  airfoils,  inasmuch  as  the  viscous  effects  reduce  the  overall  level  of  recompression  and 
bring  one  further  from  the  asymptotic  situation  of  strong  shock. 

The  non  asymptotic  two-layer  models  related  to  the  Lighthill  concepts  are  mentioned  in  the  Melnik  lecture  1 101  as  well  as  in  the  Bohning  and 
Zierep  (17|and  Inger  (181  papers.  As  in  asymptotic  theories  these  are  local  models,  usable  only  in  the  vicinity  of  shock-wave-boundary  layer 
interaction,  and  neglecting  viscosity  in  the  external  part  of  the  viscous  layer.  An  essential  difference  resides,  on  the  other  hand,  in  thedetermi 
ning  role  assigned  here  to  the  viscous  sublayer.  Recent  results  1 171(18)  are  given  in  the  form  of  parametric  studies  including  problems  of  wall 
curvature  and  incipient  separation. 


2.4  —  Computation  including  thin  layer  separations 


The  multilayer  models  available  for  turbulent  flows,  asymptotic  or  not  lose  their  validity  in  case  of  separation,  even  if  limited  to  small 
bubbles  Presently,  only  single  boundary  layer  models  make  it  possible  to  approach  the  computation  of  turbulent  separation. 


The  simplest  modeling  rests  on  the  "interactive"  solution  of  the  Prandtl  equations  and  of  a  patched  potential  flow.  It  would  be  consistent,  in 
the  laminar  case,  with  the  triple-deck  asymptotic  model.  Results  of  this  type  are  sought  for  local  interaction  phenomena,  in  incompressible 
[12]  (20)  [21 1  or  single  wave  supersonic  (27]  (28]  (29],  the  perturbation  of  the  external  flow  being  then  able  to  be  computed  in  an  explicit  way, 
or  through  Cauchy  integrals  Cebeci  et  al.  (20)  study  the  separation  at  the  leading  edge  of  an  elliptical  airfoil  at  incidence,  and  obtain  solutions 
similar  to  the  breakdown  of  small  laminar  or  transitional  bubbles.  Gleyzes  et  al.  (21J  take  as  a  basis  a  detailed  experimental  study  of  transition¬ 
al  bubbles  at  the  leading  edge  of  an  airfoil  for  defining  an  empirical  approximate  model,  allowingacomputationoftransitioninshortbubbles 
with  local  coupling  effect.  Ardonceau  et  al.  [28!  compute  the  turbulent  separation  on  a  supersonic  wedge  by  an  iterative  technique  using  an 
inverse  solution  of  the  Prandtl  equation.  The  numerical  discussion  of  the  optimal  choice  of  the  inverse  method  confirms  the  importance  of 
carrying  to  the  wall  the  coupling  condition,  a  process  known  for  eliminating  supercritical  behaviors  [1|. 


More  general  models  of  strong  interaction  with  single  boundary  layer,  capable  of  computing  separated  zones,  are  in  the  course  of  development 
for  complex  flows  around  airfoils  (1|  [9j  (11].  Transonic  computations  within  the  small  perturbation  approximation  are  mentioned  by  Wai  end 
Yoshihara  [9j,  the  shock  root  separation  remaining  however  subject  to  a  modeling  of  equivalent  viscous  wedge.  Computations  based  on  the 
solution  of  the  complete  potential  equation,  coupled  with  an  integral  method  modeling  the  discrepancy  between  the  fictitious  interactino 
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ideal  fluid  and  the  real  fluid  have  been  carried  outtLe  Balleur  (l)jThey  include  the  computation  and  positioning  of  diwymmetric  wakes  and 
of  the  separation  bubbleson  the  airfoil  or  at  the  trailing  edgei,  except  separation  bubbles  under  the  shock  in  transonic  regime.  The  progress  ob¬ 
tained  by  NAron  in  defining  a  panel  method  numerically  well  conditioned,  in  the  case  of  airfoils  with  angular  or  slender  trailing  edge,  made  it 
possible  to  adapt  the  above  method  to  the  computation  of  multicomponent  airfoils  with  multiple  separation  bubbles  in  incompressible  flow 

Ml  1- 

2.5  —  Computation  of  airfoils  and  wings-  Numerical  coupling  methods 

The  computing  methods  for  transonic  airfoils  which  have  reached  an  almost  operational  development  stage  usually  take  into 
account  the  displacement  and  wake  curvature  effects,  but  are  unable  to  predict  separations.  Elements  of  analysis  on  this  point 
are  to  be  found  in  the  Melnik  [10]  and  Le  Balleur  (11  lectures.  A  more  complete  review  has  been  given  by  the  Lock  |2] 
lecture,  presented  by  Dr  Green.  This  review  also  covers  the  much  fewer  methods  that  have  been  developed  for  transonic 
wing*  in  the  form  of  successives  iterations  between  computations  of  potential  flow  and  three-dimensional  boundary  layer.  In  view  of 
the  difficulties  still  to  be  solved,  in  this  field,  it  is  not  yet  possible  to  have  access  to  strong  coupling  methods.  However,  as 
shows  by  Firmin  [8]  on  wing-fuselage  configurations  with  the  small  perturbation  transonic  theory,  viscosity  can  already  be  taken 
into  account  to  quite  an  appreciable  degree. 

The  above  operational  computing  methods,  applicable  to  transonic  airfoils  and  wings  without  separation,  rest,  for  achieving  the 
numerical  coupling  between  ideal  fluid  and  viscous  layers,  on  a  direct,  iterative  computation  of  the  displacement  effect,  the  con¬ 
vergence  of  which  is  all  the  less  assured  as  separation  is  closer  or  the  computing  mesh  pattern  is  finer.  In  these  methods,  the  computa 
tion  stability  Is  obtained  at  best  by  relatively  empirical  techniquesof  under  relaxation.  Recent  progress,  on  the  contrary,  induce 
one  to  consider  as  a  central  element  of  strong  coupling  methods  the  mathematical  definition  of  numerically  stable  coupling  algo¬ 
rithms.  An  analysis  of  this  question,  as  well  as  a  review  of  the  main  numerical  methods  already  developed  for  strong  coupling, 
are  given  in  the  introductory  general  lecture,  Le  Balleur  [1],  In  the  steady,  two-dimensional  case,  the  direct  relaxation  methods 
can  be  stabilized  without  empiricism  by  computing  a  local  relaxation  coefficient,  while  semi-inverse  or  semi-implicit  relaxation  techniques  makes 
it  possible  to  approach  the  strong  coupling  of  separated  boundary  layers.  The  papers  of  Wai  and  Yoshihara  [9]  Le  Balleur  and 
N4ron  (11),  Veldman  (12),  Ardonceau  et  al.  (28]  complete  the  analysis  of  the  problem. 

2.6  -  Other  types  of  approach 

Apart  from  the  other  computing  methods  mentioned  in  this  Symposium,  the  paper  by  Ashurst  et  al.  |24]  provides  an 

attempt  of  evaluation  of  Lagrangian  methods  using  discrete  vortex  particles  for  simulating  numerically  the  large  scale  turbulent 
separations.  Applied  to  reattachment  behind  a  backward  facing  step,  and  restricted  to  two-dimensional  computing  assumptions, 
although  an  attractive  qualitative  simulation  is  obtained,  the  technique  does  not  seem  likely  to  match  experiments  for  deter¬ 
mining  the  Reynolds  stresses,  or  more  simply  for  reproducing  the  variation  of  reattachment  length  as  a  function  of  Reynolds 
number. 

It  would  seem  interesting,  as  may  be  suggested  by  the  Kline  general  lecture  [22] ,  to  place  the  proceedings  of  the  present  Sym¬ 
posium  in  the  wider  framework  of  the  computation  of  complex  turbulent  flows,  which  has  been  chosen  as  the  theme  of  the 
1980-81  Stanford  Conference.  This  intends  to  determine  a  range  of  detailed  experiments  covering  turbulent  flows  as  varied  as 
possible,  and  to  use  them  to  test  the  numerical  methods  presently  available. 

3.  PARTICULAR  COMMENTS 

More  specific  remarks  may  be  outlined  from  this  Symposium,  based  on  the  example  of  the  computation  of  wing  airfoils,  v  hich 
constituted,  to  a  great  extent,  the  center  of  interest  of  the  papers.  This  example  indeed  gathers  most  of  the  viscous  n  ter  ac¬ 
tion  problems  of  two-dimensional,  transonic,  external  aerodynamics,  and  constitutes  a  prelude  to  developments  of  three-dimensio¬ 
nal  computations. 

3.1  -  Use  of  the  potential  equation 

All  papers  are  based,  for  determining  the  external  ideal  fluid  flow,  on  the  solution  of  the  potential  equation,  possibly  simplified 
within  the  framework  of  the  small  transonic  perturbation  approximation.  Therefore,  progress  is  still  possible  as  regards  computation  of 
the  two-dimensional  ideal  fluid  flow.  This  remark  has  already  been  confirmed  by  a  number  of  attempts  presented  in  this  Symposium  to 
improve  the  numerical  conditioning  of  the  panel  methods  used  for  multiple  airfoils  at  low  speeds.  The  remark,  also 
applies  essentially  to  the  absence  of  use  of  numerical  methods  solving  the  Euler  equations,  probably  because  of  their  higher  cost, 
although  the  local  Mach  numbers  involved  in  transonic  viscous  fluid  around  the  incipient  separation  under  the  shock  makes 

questionable  the  potential  approximation.  This  first  limitation,  of  inviscid  origin,  should  be  kept  in  mind  in  comparisons  with 
experiments.  Furthermore,  it  has  given  rise  to  several  schools  of  thought  for  the  practical  utilization  of  transonic  potential  nume 
rical  methods,  based  on  the  exploitation  of  numerical  techniques  sometimes  conservative,  sometimes  non-conservative,  sometimes 
partially  conservative  as  regards  mass  fluxes.  This  confusion  is  all  the  more  delicate  as  it  usually  interferes  with  the  adoption  of 
a  more  or  less  elaborate  treatment  for  the  coupling  of  viscous  effects  at  the  shock  root,  the  overall  result  being  then  dependent 
on  some  compensation  of  errors  between  the  two  effects.  This  theoretical  obstacle  applies  especially  to  non-conservative  techni 
ques,  best  known  for  their  results  closer  to  experiment  in  the  absence  of  any  viscous  effect  The  conservative  techniques,  proper 
ly  coupled  with  viscous  effects  under  the  shock,  seem  nowadays  to  offer  the  safest  way  for  the  difficult  evaluation  of  drag 
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3.2  —  Use  of  integral  equations 

A  second  remark  brought  to  light  by  airfoil  and  wing  computation  is  the  general  use  of  integral  methods,  most  often  of  entrain 
ment  methods,  to  compute  turbulent  viscous  layers.  While  this  remark  does  not  apply  to  all  analyses,  especially  to  some  studies 
of  local  viscous  phenomena,  and  while  it  may  be  applied  first  with  a  view  to  reducing  the  cost  of  numerical  computation,  espe¬ 
cially  in  the  prospect  of  an  extention  to  three-dimensional  problem*,  it  is  clear  however,  that  this  choice  is  justified  by  the  qua¬ 
lity  of  the  practical  results  which  can  be  obtained  for  mean  velocity  profiles  and  for  displacement.  This  opportunity  of  numerical 
simplification  could  also  possibly  lead  to  the  computation  of  separation  bubbles,  as  well  as  to  the  incorporation  of  a  computation 
of  hystirisis  effects  on  the  turbulent  stresses  for  improving  the  modeling  of  entrainment  Moreover,  the  progress  resulting  from 
a  defect  formulation  of  the  viscous  equations  relative  to  a  superposed  fictitious  flow  of  ideal  fluid  1 1 1  emphasizes  the  fact 
that  integral  methods  actually  constitute  a  computing  model  less  restrictive  than  the  Prandtl  equations.  Their  utilization,  couplet, 
with  the  solution  of  the  fictitious  flow  of  ideal  fluid  extending  the  external  flow,  leads  indeed  to  an  approximate  representation 
of  normal  pressure  gradients  closer  to  that  given  by  the  Navier -Stokes  thin  layers  equations.  It  tolerates  in  particular  the  appea 
ranees  of  shock-waves  within  the  boundary  layers. 

3.3  -  Wake  curvature  effect 

The  influence  of  the  wake  displacement  effect  on  the  computation  of  the  global  flow  is  unanimously  acknowledged.  However, 

despite  the  apparent  consensus  on  the  importance  of  the  wake  curvature  effect,  various  approaches  have  been  proposed  according  to 
the  technique  adopted  for  defining  a  mean  curvature  of  the  streamlines.  Whereas  the  various  formulations  are  almost  equivalent 
for  the  far  wake  which  exerts  only  a  small  influence,  a  dissymmetric  treatment  of  the  near  wake  is,  on  the  other  hand,  probably 
necessary  for  correctly  performing  a  strong  interaction  computation.  This  computation  would  also  require  the  use  of  a  stable  cou¬ 
pling  algorithm,  free  of  smoothing.  It  is  probably  because  of  these  difficulties  that  a  true  consensus  has  not  yet  been  reached  on 
how  to  take  into  account  the  viscous  effects  at  the  trailing  edge  in  a  turbulent  regime,  except  in  the  limiting  case  of  an  infinite 
Reynolds  number  The  consequences  on  the  practical  evaluation  of  drag  seem  to  be  potentially  important. 

3.4  -  Shock-boundary  layer  interaction 

Opinions  also  diverge  on  the  choice  of  a  practical  treatment  of  the  shock  wave-boundary  layer  interaction  adapted  to  the  compu¬ 
tation  of  transonic  airfoils  and  wings.  To  the  already  mentioned  uncertainties  of  the  numerical  technique  adopted  for  the  external 
flow  is  added  the  difficulty  of  a  discretization  mesh  pattern  usually  too  coarse  for  this  phenomenon,  the  total  extent  of  which  is 
often  smaller  than  the  computing  mesh  size  Although  it  is  generally  admitted  that  the  overall  consequences  of  this  phenomenon  on 
the  flow  are  less  determining  than  those  resulting  from  the  trailing  edge  region,  at  least  two  requirements  must  be  fulfilled  by  the 
computation  :  the  approximate  viscous  computation  should,  on  the  one  hand,  modify  fundamentally  the  evolution  of  the  recom¬ 
pression  under  the  shock  and  reduce  its  level.  This  precludes  the  direct  use  of  an  asymptotic  model  with  normal  shock  and  may 
require  in  some  cases  the  addition  of  a  technique  of  artificial  viscous  wedge  ;  the  approximate  computation  should,  on  the  other 
hand,  estimate  the  overall  evolution  of  the  boundary  layer  with  a  precision  sufficient  for  validly  computing  the  trailing  edge  re¬ 
gion.  The  future  development  of  solutions  on  a  finer  scale,  guaranteeing  the  computation  and  positioning  of  the  phenomenon,  seems 
however  indispensable. especially  ir.  case  of  a  local  separation. 

3.5  —  Operational  codes  and  the  necessity  of  experimental  checking 

In  spite  of  these  unsolved  difficulties,  transonic  airfoil  computing  methods  based  on  the  solution  of  the  potential  equation  cou¬ 
pled  with  integral  methods  for  the  viscous  layers  have  reached  an  operational  level  of  practical  prediction,  at  least  in  the  absence 
of  separation.  The  performance  obtained  is  better  than  that  of  the  methods  based  on  direct  solutions  of  averaged  Navier -Stokes 
equations,  as  regards  cost,  quality  and  mesh  size.  These  encouraging  results  in  steady  two-dimensional  flow  should  neither  hide  the  above 
mentioned  insufficiencies,  nor  the  difficulty  of  objective  comparisons  between  computations  and  transonic  experiments  which  appear 
at  this  level  of  approximation.  It  would  be  particularly  desirable  to  have  available  reference  experiments  strictly  two-dimensional,  (i.e. 
almost  exempt  of  wall  interference,)  comprizing  a  very  accurate  description  of  the  boundary  layers  behind  the  transition 
onset.  Thus,  the  comparison  to  the  computation  could  be  carried  out  for  a  given  Mach  number  and  angle  of  attack,  rather  than 
by  the  very  questionable  bias  of  an  adjustment  of  the  overall  lift. 

3.6  -  Outlook  for  running  and  further  works 

Recent  progress  achieved  on  the  mathematical  models  coupling  ideal  fluid  and  viscous  layers,  as  well  as  on  the  numerical  al¬ 
gorithms  making  it  possible  to  solve  the  corresponding  strong  coupling  problems,  has  reached  a  stage  allowing  the  application 
of  the  above  mentioned  two-dimensional  methods  to  regimes  of  thin  layer  separation.  In  this  domain,  which  corresponds  to 
high  lift  coefficients,  multicomponent  airfoils  and  transonic  buffeting,  an  increased  sensitivity  to  turbulent  modelings  as  well 
as  to  the  treatment  of  the  transition  zone,  should  however  be  expected  in  the  practical  performance  of  the  computing  me¬ 
thods.  Furthermore,  a  theoretical  link  between  these  models  of  separation  and  reattachment  calculation,  and  the  simpler  mo 
dels  with  isobaric  jet  lines  of  extensive  separations,  should  be  sought  in  the  future. 

Although  the  essential  part  of  the  study  remains  to  be  carried  out,  the  extention  of  strong  coupling  methods  based  on  a  solution 
in  ideal  fluid  of  the  external  flow  and  on  integral  methods  for  the  viscous  layers  should  be  ab'e  to  make  possible  the  approach 
of  the  practical  computation  of  three-dimensional  or  unsteady  flows  with  thin  layer  separations.  To  this  end,  the  newest  theore- 
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tics)  and  numerical  obstacle  to  be  overcome  probably  consists  of  imagining  a  viscous  model  of  strong  interaction  describing  the 
shedding  and  rolling  up  of  the  wake  vortex  sheets  issued  from  slender  bodies. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  Symposium  clearly  proves  that  the  effort  devoted  these  last  years  to  numerical  methods  of  ideal  fluid  computations  may  be 
highly  enhanced,  from  a  practical  viewpoint,  by  the  addition  of  a  complementary  interactive  solution  of  the  viscous  layers. 

The  work  is  well  advanced  as  regards  steady  two-dimensional  flows,  using  a  potential  approximation  of  the  exterral 
flow.  The  state  of  the  art  in  this  field  lies  at  an  intermediate  level  between  weak  coupling  methods,  isoberic 
modeling  methods  for  extensive  separations  and  strong  coupling  methods,  these  last  ones  being  capable  of  representing  the  com¬ 
plete  set  of  real  phenomena  of  viscous  interaction,  at  least  in  a  qualitative  fashion. 

The  latest  progress  tends  toward  strong  coupling  methods,  with  the  joint  development  of  mathematical  models  coupling  ideal 
fluid  and  viscous  layers,  as  vwll  as  of  the  numerical  algorithms  necessary  for  their  exact  coupled  solution.  Although  uncertain¬ 
ties  remain  concerning  problems  of  trailing  edge  or  shock  wave-boundary  layer  interaction,  operational,  high  performance  me¬ 
thods  exist  for  transonic  airfoils  ;  one  may  observe  in  particular  the  appearance  of  new  methods,  capable  of  computing  sepa¬ 
rations  in  thin  layers,  or  trailing  edge  separations. 

The  implementation  of  the  solutions  of  the  complete  Euler  equations  for  the  external  flow  might  seem  desirable,  for  elimina- 
ting  on  the  one  hand,  the  difficult  choice  between  conservative  and  non-conservativc  potential  techniques,  and  on  the  other  hand 
for  extending  the  field  of  application  of  the  strong  coupling  methods  to  problems  of  transonic  internal  aerodynamics. 

The  field  of  unsteady  or  three-dimensional  flows  is  far  less  developed.  However,  the  first  results  obtained  by  a  solution  of 
the  ideal  fluid  flow  by  the  potential  approximation,  and  by  addition  of  a  coupled  computation  of  the  viscous  layers  by 
means  of  integral  methods,  seem  to  prove  the  possibility  of  generalizing  the  two-dimensional  models. 

The  particularly  encouraging  performance  of  the  strong  coupling  methods  for  practical  applications,  for  which  they  consti¬ 
tute  a  tool  of  approximate  numerical  simulation  of  the  Navier-Stokes  equations  for  thin  layers,  leads  one  to  recommend  for 
the  future  associating  more  systematically  and  more  closely  the  development  of  these  numerical  methods  with  techniques  for 
solving  the  ideal  fluid  equations. 

In  the  same  manner,  future  and  complementary  developments  of  strong  coupling  methods,  as  well  as  global  techniques  for 
solving  the  Navier-Stokes  equations,  should  probably  be  more  closely  associated.  This  is  the  case,  for  problems  where  the 
"Navier-Stokes"  computing  domains  are  limited  to  viscous  regions  only.  Moreover,  it  appears  possible  to  develop  new  algorithms 
for  the  Navier-Stokes  equations,  based  on  numerical  techniques  derived  from  strong  coupling  methods,  at  least  within  the  thin 
layer  approximation. 
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S  -2  § 

*"*»  o  is) 

■a  2  5 

■S  2  S 

>  «5  c 
•S 

s  I  § 


O  C  ^ 

°  ff 

5  "o  v,  * 
o  4J  a> 

S  </l  »— 
eg 

e  X>  J4  , 

I  </>  H  - 

Ilf 

a  2  c 
«  s  o  • 
"  a  °  , 
.2  <u  o 
—  ^  .c 


§52 

MS 

**  «* 

<c* 

(/)  a.  -C 


u  •  — 

x:  r- 
—  r  i  -a 
^ 

eg  <*  *c 

.<2 

©  3  35 
00—3 
^  Sj  i. 


^  — 

5  .2  ! 

a.  u  S 
■2  2  g 


x  c  r  ’ 

Ilf 

C.  2  fc 

="  B  9  ■ 


CJ  Ct  Q 

2  £  o£ 

t2  "o  < 

uS'j 

<  c  S 

on  a.  ~ 
T  o-  c 


fl  tf)  -C 

^  c/s 

0  3  2 
00—3 
O'  ,J  ? 

zL  a>  c* 


5  2  3 
0  «—  <*• 
c  ~  .ts 
-.SB 

eg  G-. 

§  "S  « 

.X  O  X) 


.§13  | 

i 

—  t/3 

C/3  c/3 

V  OJ  »— « 

■CT  S'  <u 

m  5?  2 

Soo 

Si  ft.  £ 


o  c  ° 
o  8 

0—3 
•=  2  o 

3  •£ 
jz  c.  2 

•s  e  e 
s  3  £ 

M3 

o  *c  u— 
c:  w  • - 
C  3 
«  ’  “  ^ 
§  i  «> 

.2  O  X3 

V»  c 

5  ?i 

•i  1  | 

•V  g  * 

—  C/3 

C/3  ^ 

5  p  o 


iL  ~  c  : 

«u  •->  o  Yi 
—  c  ■—  ■*■> 

-|So 

tcS” 

£  s  ,s  .£ 

c  > 
•o  ?  C  Z 

V3  —  O  2 
01 


Academy  des  Etats-Unis,  a  Colorado-Springs,  Colorado,  USA.  Lcs  4  confdrences,  21 

communications  et  la  Discussion  de  Table  Ronde  presentees  lors  du  Symposium  sont  This  Advisory  Report  was  produced  at  the  request  of  the  Fluid  Dynamics  Panel  of 
publiees  dans  i’AGARD  Conference  Proceedings  CP-29 1 ,  date  de  f£vricr  1981 .  AGARD. 


